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We provide a brief overview of some current experimental and theoretical issues of 
(-H ' heavy quark production. 
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1 Introduction 



The production of heavy quarks in high energy processes has become an increasingly 
^. ' important subject of study both theoretically and experimentally. The theory of 

l/*~) , heavy quark production in perturbative Quantum Chromodynamics (pQCD) is 

UN ' more challenging than that of light parton (jet) production because of the new 

^j| , physics issues brought about by the additional heavy quark mass scale. The correct 

(—} ' theory must properly take into account the changing role of the heavy quark over 

Qs, . the full kinematic range of the relevant process from the threshold region (where 

Q\ ' the quark behaves like a typical "heavy particle" ) to the asymptotic region (where 

the same quark behaves effectively like a massless parton). With steadily improving 

O ' experimental data on a variety of processes sensitive to the contribution of heavy 

' , quarks (including the direct measurement of heavy flavor production), this is a very 

rich field for studying different aspects of the QCD theory including the problems of 

C^ , multiple scales, summation of large logarithms, subtleties of renormalization, and 

higher order corrections. We shall briefly review a limited subset of these issues.^] 

X 

2 The Factorization Theorem 

Perturbative calculations for heavy quark production are performed in the context 
of the factorization theorem expressed below in the commonly used form: 

o a ^ c = f a ^ h {x, m 2 ) ® a b ^ c (Q 2 /n 2 , M 2 H /v 2 ,a s (p)) + 0(A 2 QCD /Q 2 ) (1) 

While the factorization was originally proven for massless quarksp the theorem has 
recently been extended by Collinsa to incorporate quarks of any mass, including 
"heavy quarks." (Note, we have explicitly retained the Mfj dependence in a.) It 
is important to note that the corrections to the factorization are only of order 
Aq CD /Q 2 , and not Mfj/Q 2 , even {or the case of general quark masses. 

The factorization theorem can also be expressed as a composition of i-channel 
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Figure 1. Basic processes for DIS heavy quark production, a) Q(oq) flavor excitation: 7 + c — » c; 
b) O(aJ) flavor creation: 7 + g — * c + c; c) 0(a\) flavor excitation: 7 + c — > c + g; d) 0(a\) 
light-quark (q) fragmentation: (7 + q — » g + g) (g) (g — > c) . 



two particle irreducible (2PI) amplitudes:^ 
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Here, C represents the graph for a hard scattering, K represents the graph for a 
rung, T represents the graph that couples to the target, and Z represents a collinear 
projection operator. The first term in brackets roughly corresponds to the hard 
scattering coefficient function a, and the second term to the parton distribution 
function (PDF), /. Note that these two terms only communicate through a collinear 
projection operator, Z. Part of the effort in generalizing the factorization theorem 
for the case of massive quarks involves constructing the proper Z, and demonstrating 
that terms containing (1-Z) are power suppressed. However, once Z is determined, 
Eq.ph yields an all-orders prescription for computing both the hard scattering 
coefficient (a) and the parton distribution function (/). A calculation using this 
formalism was first performed by ACOTq for the case of heavy quark production 
in deeply inelastic scattering, and we now examine this process in detail. 



3 Heavy Quark Production in DIS 

Several experimental groupsB have studied the semi-inclusive deeply inelastic scat- 
tering (DIS) process for heavy-quark production t\ + N — > £2 + Q + X. New data 
from HERA investigates the DIS process in a very different kinematic range from 
that available at fixed-target experiments. This perception has changed the way 
that we compute the semi-inclusive DIS heavy quark production. Traditionally, the 
heavy quark mass was treated as a large scale, and the number of active parton 
flavors was fixed to be the number of quarks lighter than the heavy quark. In this 
scheme, the perturbation expansion begins with the 0(a l s ) heavy quark creation 
fusion process jg — > cc, (c/., Fig. fib). We refer to this approach as the Fixed Flavor 
Number (FFN) scheme since the number of flavors coming from parton distributions 
is fixed at three for charm production]] 



b l must necessarily leave out many details here; for a precise treatment, see Collinstj. 

c The necessary diagrams have been computed to 0(a 2 s ) by Smith, van Neerven, and collaborators, 

cf., Ref. B 
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Figure 2. The scaled cross section for DIS heavy quark production as a function of the quark mass 
m H . 



More recently, a Variable Flavor Number (VFN) scheme (ACOTcl) has been 
proposed which includes the heavy quark as an active parton flavor with non-zero 
heavy quark mass. In this case, the perturbation expansion begins with the 0(oq) 
heavy quark excitation process jc — > c, (c/., Fig. |l|a). The key advantages of this 
scheme arefl 

1. By incorporating the heavy quark into the parton framework, the composite 
scheme yields a result which is valid from threshold to asymptotic energies; in 
contrast, the FFN scheme contains unsubtracted mass singularities which will 
vitiate the perturbation expansion in the m c — > or E — > 00 limit. 

2. Because the composite scheme resums the large logarithms appearing in the 
FFN scheme into the parton distribution functions, it includes the numerically 
dominant terms of the 0(a 2 s ) FFN scheme calculation in a 0(a\) calculation. 

In effect, the VFN scheme subsumes the FFN scheme. To illustrate this fact with a 
concrete calculation, in Fig. pi we plot the cross section for "heavy" quark produc- 
tion as a function of the quark mass.fi This figure clearly shows the three important 
kinematic regions. 1) In the massless region, where ran <C Q, the ACOT VFN re- 
sult reduces precisely to the massless MS result. 2) In the decoupling region, where 
itih 3> Q, this "heavy quark" decouples and its contribution vanishes. 3) In the 
transition region, where run ~ Q, this (not-so) "heavy quark" plays an important 
dynamic role. While the FFN scheme is appropriate only when ma ~ Q, we see 
that the VFN scheme is valid throughout the full kinematic range.f] 

This point is also illustrated in a calculation by KretzerQ (c/., Fig. y) which 
shows the partial contributions to the charged current i^ ™" v\ In this figure, 



d To be specific, we have computed single quark production for a photon exchange with x = 0.1, 

(j, = Q = 10 GeV, and the cross section is in arbitrary units. 

e Buza et ah, have determined the asymptotic form of the heavy quark coefficient functions which 

are then used rto determine the threshold matching conditions between the three- and four-flavor 

shemes, Ref. Q.i-iThorne and Roberts have a similar scheme with slightly different matching 

conditions, Ref. U- 

^Kretzer and Schienbein have performed the first calculation of the 0(0:5) quark initiated process 
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Figure 3. The contributions to DIS charged current inclusive F^ 



vs. /^. For each separate 



contribution, the thick lines are the MS rpsult (m s 
with m a = 0.5 GeV. From Kretzer, Ref. 13. 



0), and the thin lines are the ACOT result 



each line is actually a pair of lines: the thin lines represents the result for F£ harm 
using the ACOT scheme with m s — 0.5 GeV, and the thick lines regularize the 
strange quark with the massless MS prescription. (The charm mass is, of course, 
retained.) The fact that these two calculations match so closely (particularly in 
comparison to the /z-variation) indicates: 1) the ACOT scheme smoothly reduces 
to the desired massless MS limit as run — » 0, and 2) for ttih ^ A-qcd we see that 
the quark mass no longer plays a dynamic role in the process and becomes purely 
a regulator. 



4 Heavy Quarks and Extraction of s(x) 

A topic closely related to DIS charm production is the extraction of the strange 
quark distribution.liiP In principle, we can extract s(x) by comparing DIS neutral 
and charged current data. To leading order, we have: 



F 2 NC 

f£° 



_5_ f 3 (s + s)- (c + c) + 
18 1 ~5 VTd 



(3) 



While the individual Fi structure functions are measured precisely (c/., Fig. 0),M 
this approach is indirect in the sense that small uncertainties in the larger valence 
distributions will magnify the uncertainty on the extracted s{x). 

A direct method of obtaining s(x) is to use the neutrino induced di-muon pro- 
cess: VpN — > (i~cX with the subsequent decay c — » sy^v^. Here, the di-muon 
signal is directly related to the charm production rate, which goes via the process 
W + s — > c at leading order. The method has the advantage that the signal from 
the s-quark is not a small effect beneath the valence process. 

complete NLO experimental analysis was performed using the CCFR data 



set.trl The recently collected data from the NuTeV experiment will provide an 



for general masses and general couplings, Ref. Q. i— i 

9 For a comprehensive review, see Conrad, Shaevitz, and Bolton, Ref. E3. 




Figure 4. Comparison of F2 from charged and neutral current DIS. From Seligman, et ah, Ref. 



B. 



opportunity to extend the precision of these investigations still further. Ell Their 
high intensity sign-selected neutrino beam and the new calibration beam allows 
for lame improvement in the systematic uncertainty while minimizing statistical 



errors. I 



5 Hadroprodution of Heavy Quarks 



We now tups to the hadroproduction of heavy quarks, and discuss how the method 
of ACOTqO is used to provide a dynamic role for the heavy quark parton. We con- 
centrate mostly on b-productioi] af fh p TWatmn for definiteness, and present typical 
results for b quark productionJdOllallj (See the paper by D. Fein, this meeting.Ej) 
Fig. |k shows the scaled differential cross section vs. pr for b production at 1800 
GeV for the leading order (LO) calculations. The heavy creation (HC) processn 
(gg — > bb) represents the LO contribution to the fixed-flavor-number (FFN) scheme 
result. The heavy excitation (HE) process (gb — > gb) plus the HC term represents 
the LO contribution to the variable-flavor-number (VFN) scheme result. The pair 
of lines for each result shows the effect of varying fi. In a similar manner, Fig. |3|b 
shows the total FFN and VFN results]] 

Two interesting features are worth noting. 1) Examining Fig. gja we observe 
the HE contribution is comparable to the HC one, in spite of the smaller &-quark 
PDF compared to the gluon distribution. Closer examination reveals that two 



h In this section we let g represent both gluons and light quarks, where applicable. Therefore, the 
HC process described as gg — > bb also includes qq —* bb. 

8 The formifia,ble calculations of the NLO gg —* bb rerecess were computed by Nason, Dawson, and 
Ellis (Ref. EU), and also by Beenakker et al., (Ref. E3). These calculations were implemented in a 
Monte Carlo framework (including correlations) by Mangano, Nason, and Ridolfi , (Ref. Ej). 
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Figure 5. The scaled cross section (nb-GeV 3 ) vs. pt for the a) LO-FFN and LO-VFN con- 
tributions, and b) TOT-FFN and TOT-VFN contributions. For each contribution, we choose 

fj, = [Mt/2,2Mt], with Mt = \/' m fj + Pj" to § au § e the //-variation. From Ref. £3. 



effects contribute to this unexpected result: a larger color factor and the presence 
of t-channel gluon exchange diagrams for the HE process, as compared to the HC 
process. 2) The LO-VFN (=HC+HE) contributions (Fig. ||a) (tree processes) give 
a reasonable approximation to the full cross section TOT-VFN (Fig. gb); thus, the 
NLO-VFN correction is relatively small. This is in sharp contrast to the familiar 
FFN scheme where the TOT-FFN term is more than twice as large as the LO-FFN 
(=HC). This is, of course, an encouraging result, suggesting that the VFN scheme 
heavy quark parton picture represents an efficient way to organize the perturbative 
QCD series. 

In Fig. ga, we also observe that while the TOT-VFN result provides minimal 
/^-variation for low pr, the improvement is decreased for large px- This may be, in 
part, due to that fact that the TOT-VFN result shown here is missing the NLO-HE 
process gb — ► ggb since this calculation, with masses retained, does not exist. In 
a separate effort, Cacciari and GrecoEJ have used a NLO fragmentation formalism 
to resum the heavy quark contributions in the limit of large pt- This calculation 
effectively includes the massless limit of the gb — » ggb contribution (omitted above); 
the result is a decreased /i- variation in the large pt region. Recently, this calculation 
has been merged with the massive FFN calculation by Cacciari, Greco, and Nason, 
(Rcf. Ej); the result is a calculation which matches the FFN calculation at low pr, 
and takes advantage of the NLO fragmentation formalism in the high pt region. 

Not only do the VFN and FFN schemes yield different cross sections, but the 
differential distributions can be quite distinct. Specifically, in Fig. |fj|a we display 
the rapidity distribution for TOT-VFN as compared with the TOT-FFN result. 
We observe that the VFN scheme yields a broader rapidity distribution than the 
FFN scheme; in part, this is expected as the VFN result includes a t-channel gluon 
exchange process (a + b — ► g + b) which can give an enhanced contribution in the 
forward direction.[] It is interesting to compare this result with the data; Fig. gb. 
shows the comparison with the TOT-FFN result. In the central region, the TOT- 



■'Note, the increased cross section in the forward region is not guaranteed a prioi. Only after the 
collincar singularity has been subtracted (as we have done for Fig. Ha) can the sign of the effect 
be determined. 
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Figure 6. a) Comparison of the total cross section da/dp^/dy in the VFN and FFN formalism 
vs. y for b production at 1800 GeV with pt = 20 GeV and (j, = Mt- The curves are scaled rby 
the total FFN cross section to facilitate comparison of the relative magnitude. From Ref. lij. 
b) The measured differential muon cross sections from b production and decay as a function of 
y. The Solid curve is a TOT-FFN predict*©* with dashed curves representing the theoretical 
uncertainties. Preliminary results from Ref. I 18 l . 



FFN theory is a factor of ~2 below the data; this increases to about a factor of ^3 
in the forward region. The TOT- VFN scheme increases (as compared to the TOT- 
FFN result) at large rapidity, as do the data. This observation suggests that the 
VFN scheme (which includes the flavor excitation process g + b — ► g + b) , provides 
a mechanism that can help resolve the y shape discrepancy. 

6 Massive vs. Massless Evolution 

In a consistently formulated pQCD framework incorporating non-zero mass heavy 
quark partons, there is still the freedom to define parton distributions obeying either 
mass- independent or mass-dependent evolution equations. With properly matched 
hard cross-sections, different choices merely correspond to different factorization 
schemes, and they yield the same physical cross-sections. We demonstrate this 
principle in a concrete order a s calculation of the DIS charm structure function.Ea 
In Fig. [7] we display the separate contributions to F^ harm for both mass- independent 
and mass-dependent evolution. The matching properties are best examined by com- 
paring the (scheme-dependent) heavy excitation F^ E and the subtraction Fjf 75 
contributions of Fig. |^a. 

We observe the following. 1) Within each scheme, F% E and F^ 173 are well 
matched near threshold, c/., Fig. [7|a. Above threshold, they begin to diverge, but 
the difference (F 3E — F 2 SC/S ), which contributes to F^ 07 ', is insensitive to the 
different schemes. 2) It is precisely this matching of F^ E and F^ UB which ensures 
the scheme dependence of F™ T is properly of higher-order in a s , (c/., Fig. Mb). 

This matching is not accidental, but simply a result of using a consistent renor- 
malization scheme for both F^ E and F^ . To understand this we expand these 
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Figure 7. F 2 vs. ft for DIS c-production. a) Ff fc , F£ UB and the difference (Kf fc - F^ ua ). 
The solid curves are for the mass-independent evolution scheme, and the dashed curves are for 
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the mass-dependent evolution scheme, b) F^ OT and F^ E 

the mass-independent evolution and mass-dependori*l| evolution for F"^ OT is higher order and 

comparable or less than the /^-variation. From Ref. Ej. 



terms near threshold (fi ~ ran) where the rnn/Q terms are relevant: 
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Here, the prescript i? specifies the renormalization scheme. From these relations, it 
is evident that i r ^' B and F.f 7 - 8 will match to 0(a^) so long as a consistent choice 



or renormalization scheme R is made for the splitting kernels, ± g 



(i) 



R P^c, This is 

the key mechanism that compensates the different effects of the mass-independent 
vs. mass-dependent evolution, and yields a otot which is identical up to higher- 
order terms. The lesson is clear: the choice of a mass-independent MS or a mass- 
dependent (non-MS) evolution is purely a choice of scheme, and becomes simply 
a matter of convenience- there is no physically new information gained from the 
mass- dependent evolution. 



7 Conclusions 



We have provided a brief overview of some current experimental and theoretical 
issues of heavy quark production. The wealth of recent heavy quark production data 
from both fixed-target and collider experiments will allow us to to extract precise 
measurements of structure functions which can provide important constraints on 
searches for new physics at the highest energy scales. As an important physical 
process involving the interplay of several large scales, heavy quark production poses 
a significant challenge for further development of QCD theory. 

We thank J.C. Collins, R.J. Scalise, and W.-K. Tung for valuable discussions, 
and the Fermilab Theory Group for their kind hospitality during the period in 
which part of this research was carried out. This work is supported by the U.S. 
Department of Energy and the Lightner-Sams Foundation. 
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